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The cos 20 azimuthal asymmetry of unpolarized dilepton production at the Z-pole 
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We calculate the Boer-Mulders effect contribution to the cos 2(f) azimuthal asymmetry of unpolarized dilepton 
production near the Z-pole. Based on the tree-level expression in the transverse momentum dependent factor- 
ization framework, we show that the corresponding asymmetry near the Z-pole is negative, which is opposite 
to the asymmetry in the low Q^ region, dominated by the production via a virtual photon. We calculate the 
asymmetry generated by the Boer-Mulders effect near the Z-pole at RHIC, with ^fs = 500 GeV. We find that the 
magnitude of the asymmetry is several percent, and therefore it is measurable. The experimental confirmation 
of this sign change of the asymmetry from the low Q^ region to the Z-pole provides direct evidence of the chiral 
odd structure of quarks inside an unpolarized nucleon. 
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I. INTRODUCTION 

The cos 2ip angular distribution of dilepton production in 
unpolarized hadron collisions hih2 — > (^{^ + X belongs to 
the remaining challenges which need to be understood from 
QCD dynamics [1]. According to the Lorentz structure of 
the hadronic tensor, one can write down the differential cross- 
section of dilepton production as [2, 3]: 



dcr 



= Wr(l + cos^ e) + Wl{1- cos^ 0) 
h-Wa sin 26* cos + Waa sin^flcos^ 



(1) 



Here q is the virtual photon or Z boson's four momentum, and 
dQ. = d cos 0d(p is the solid angle of the lepton { in terms of 
its polar and azimuthal angles in the center-of-mass system 
(c.m.s.) of the lepton pair. The coefficient functions Wt,l,a,aa 
depend on the invariant mass Q, transverse momentum qr, 
and the rapidity y of y* /Z. After the solid angle dQ. is inte- 
grated over, the diff'erential cross-section with respect to q has 
the form 



^ = faWr^W,). 



(2) 



Therefore, the angular distribution of the lepton pair is defined 

as 



dN do- I dcr 

dQ. ~ dQd^qj '¥q ' 



(3) 



Equivalently, another often used convention for the dilepton 
angular distribution is 



dN_ 

dn 



1 



471^-1-3 



I + A cos 6 + fi sin 2 cos c 



+ - sin 6* COS. 
2 



Comparing (1) and (4) yields following relations 

Wt - Wl Wa 21Vaa 



A^ 



Wt + Wl 



/^ = 



Wt + Wl 



Wt + Wl 



(4) 



(5) 



Of particular interest are the angular distribution given by 
the A and v terms. To Os order of perturb ative QCD, a calcu- 
lation [4] in collinear factorization showed that these coeffi- 
cients satisfy: 



2// + /I - 1 = , 



(6) 



the so-called Lam-Tung relation [3, 5], which has attracted 
considerable attention. Fixed-order pQCD calculations [6] at 
order a^, as well as QCD resummation calculations [7, 8] 
to all orders in collinear factorization, indicate that viola- 
tions of (6) are very small. However, early measurements 
on K-N -^ y* +X ^ rr +X processes by the NAIO [9] 
and El 65 [10] Collaborations at ^/s ^ 19 and 23 GeV, show 
large positive values of v, near 30%, indicating a sizable vi- 
olation of the Lam-Tung relation. The relation has also been 
tested in pp and pd Drell-Yan processes by the E866/NuSea 
collaboration [11, 12] at V^ = 38.7 GeV, and very recently in 
pp ^ y*/Z + X ^ tr +X by the CDF Collaboration [13] 
at V^= 1.96 TeV. 

Several attempts have been made to interpret these data, in- 
cluding QCD vacuum effects [14, 15] and higher-twist mech- 
anisms [16, 17]. In Ref. [18] Boer demonstrated that the 
product of two transverse momentum dependent (TMD) Boer- 
Mulders functions hf{x,p^) [19] can produce unsuppressed 
cos2(/> asymmetries that correspond to a violation of the 
Lam-Tung relation. Several theoretical and phenomenologi- 
cal studies [20-31] along this direction have been put forward. 
Those studies are mainly concentrated on the Q region much 
lower than the Z mass, where the lepton pair is produced via 
a virtual photon. In this paper, we will study the phenomenol- 
ogy of the cos 2(p asymmetry in the Z-pole region. We will 
show that the behavior of the cos 20 asymmetry coming from 
the Boer-Mulders effect in the Z mass region is very different 
from that in the low Q region. 




FIG. 1: Angular definitions of unpolarized Drell-Yan process in tire 
lepton pair center of mass frame. 



II. DESCRIPTION OF cos 20 ASYMMETRIES IN TERMS 
OF BOER-MULDERS FUNCTIONS 



The process we consider here is the dilepton production via 
a y*/Z boson in the unpolarized hadron collision: 

h{P,) + h2(P2) ^ f/Z(q) + X^ e{€) + £-((') +X. (7) 

The angular distribution coefficients /I,//, v (or W^.l.a.aa) in 
this process are generally frame dependent. In the following 
we will use the CoIIins-Soper (CS) frame [32], as shown in 
Fig. 1. This frame is the c.m.s of the dilepton, and in addi- 
tion the z-axis is chosen to be along the bisector of momenta 
Pi and -P2. In principle one can also choose the Gottfried- 
Jackson (GJ) frame [3]. The advantage of CS frame in our 
study is that in this frame the Lam-Tung relation is rather in- 
sensitive to higher-order corrections [6] and resummation ef- 
fects [7]. Besides, it is more widely used in theoretical and 
experimental studies, so that the comparison with other works 
is straight forward. 

In the TMD factorization framework, the cos 20 angular de- 
pendence arises from the coupling of the Boer-Mulders func- 
tions, which depend on the intrinsic transverse momentum, 
and appear in the decomposition of the TMD parton correla- 
tion function for an unpolarized hadron [19]: 



^{x,pt) = - 



fi(x,p\)ili+ + h\(x,p]:)- 



ypjTl_^ 



M 



(8) 



Here n+ - (0, 1,0^) is a light-like vector expressed in the 
light-cone coordinates, in which an arbitrary four-vector a is 
written as {a^ ,a* ,aT), with a^ - (a" ± a^)/y/2 and aj - 
(a\a^). The Boer-Mulders function describes the transverse 
polarization of the quark inside a unpolarized nucleon, and 
thus is chiral-odd. Despite its time-reversal odd nature, hf 
can be nonzero, due to inital/final state interactions [33-37] 
between the struck quark and the spectator of the nucleon. 
These studies have motivated the model calculations, as well 
as the lattice analysis, of the Boer-Mulders functions for nu- 
cleon [38-51] and pion [20, 21, 52-54] 

The detailed derivation of the angular dependent differential 
cross-section of reaction (7) in the TMD factorization frame- 
work has been given in Refs. [18, 55]. Here we write down 
the final expression, after taking into account both of photon 



and Z boson contributions: 

do-(hih2 -^ fix) 



iKmFl 



uu 



dQ.dxidx2Cp'qT ^Q^ 

+ [K-iiO) cos 2(p + K4(9) sin 2(p]F^J'^ 



(9) 



in which only the unpolarized production is included. In the 
above equation there are two structure functions contributing 
to the cross-section, which have the form: 



Flu = Z^[/i'/i']' 



q,q 



-,20 



Z^ 



(2h prh- kj-pT ■ kj) ' ' 
^ '^ ' M1M2 



The vector h = qr/Qr- , and we have used the notation 



/ 



ri ■ ■] = d^pTd^krdHpT + kT- qT)l ■ ■]• 



(10) 
.(11) 

(12) 



Therefore (10) and (11) represent tree-level parton model re- 
sults. 

The coefficients Kj in front of structure functions have the 
form [18] 

Km = j(l+cos^e)[el+2g'yeag''yXi+c'Ax2] 



4 
cosO 



i'^SA^aglXl + C^C^Xl] , 



Km = i sin2 [el + 2g'yeagUi + c[c^2X2] , 

Km = \^^in^ 0[2g'yeaglx^\ 
which contain the combinations of couplings 

4 = {si^-gi^), j^i or a 

4 = Ksi- 



(13) 
(14) 
(15) 

(16) 
(17) 
(18) 



The vector and axial-vector couplings to the Z boson are given 
by: 



7^3^ - 2 QJ sin^ 6^, 



SA ■'3' 



(19) 
(20) 



where Q^ denotes the charge and Ti^ the weak isospin of par- 
ticle j (for example, T^ = -1-1/2 for j - u and T^ = -1/2 for 
j - e^ , d, s). The Z-boson propagator factors are given by 



1 



Xi = 

X2 = 
X3 = 



Q'iQ'-M') 



z"^z 



sm\29w) (Q^ - M2)2 + r|M: 

1 e' 

Y] 

sm\2ew) Q^ - Ml'^ 



(21) 
(22) 
(23) 



The first term in Eq. (9) is azimuthal independent. It gives 
the (1 + cos^ 9) and cos 6 angular dependence while the later 
one vanishes after integration upon the polar angle 6. The 
structure function F^^^ therefore corresponds to Wt given in 
(1). The second term has a cos 20 azimuthal dependent term 
which contributes to the asymmetry v '. As shown in (11), it 
arises from the product of the transverse momentum depen- 
dent functions hf from each hadron. 

One important feature implied by (11) is that the cos 2(f) de- 
pendence contributed by the Boer-Mulders effect shows up in 
the pure electro-weak process without the need of QCD radia- 
tion. It will give a sizable contribution to Waa at low qr- Since 
the same effect cannot contribute to Wl, its presence violates 
the Lam-Tung relation in the low qj region. If we express the 
structure function Waa as the sum of the perturbative QCD 
effect and the Boer-Mulders effect 



WAA = W?f" + Wf^ 



'AA 



'AA 



(24) 



where W^^ is proportional to fJ"^, 
the coefficients yields 



then the combination of 



2V + A-1 



4(W,r + ^aT) 



Wt-Wl 
Wt + Wl ^ Wt + Wl 



- 1 



Wt + Wl 



= 2v 



,BM 



(25) 



The minus sign in the above equation comes from the odd 
permutation of gamma matrices in calculating the hadronic 
tensor coming from the Boer-Mulders effect: 



Tv{(TafiV^(Tp^V') 



(29) 



where V^ - gvj'^ + gAjsj'^ denotes the Z-boson-fermion ver- 
tex. 

At low Q the dilepton production via virtual photon dom- 
inates. However, as the dilepton mass approaches the Z bo- 
son mass, the production via Z-boson becomes important and 
dominates over that via virtual photon (the second term in 
(14) is suppressed in the Z mass region compared with the 
third term). Thus one expects to observe the sign reversal of 
yBM (-Qj.^ equivalently, the combination 2v + A - I that can be 
measured by experiment) from the low Q region to the Z-pole 
region. The sign reversal of v^^ is due to the chiral-odd na- 
ture of the Boer-Mulders function, and is a special feature of 
the Boer-Mulders effect as a source of the Lam-Tung relation 
violation. Therefore, the experimental detection of the sign 
reversal of 2v -H /I - 1 is a clear evidence on the existence of 
chiral-odd quarks inside the unpolarized nucleon. The exper- 
iments can be conducted at hadron colliders with unpolarized 
beams, and thus this allows to study the spin structure of the 
nucleon without polarized beams. 



Here v**^ denotes the cos 20 asymmetry contributed by the 
Boer-Mulders effect, in analogy with the definition of v in (5). 
In the above equation we have used the Lam-Tung relation 



2wff^ - Wl 







(26) 



The » sign follows the fact that higher order perturbative con- 
tributions still give a very small contribution. As explained 
previously, this contribution can be minimized by choosing 
the CS frame. At low qr, there is Wl <k Wt, then we can 
arrive at following expression for v**': 



y'^'igr, 



y,Q) 



2W 



BM 



2F] 



2(f> 



Wt 



F] 



(27) 



uu 



The approximation in the above equation comes from the tree- 
level expressions (10,1 1) for F\jjj and F^^, that is, we do not 
consider the soft factor in the TMD factorization formula. We 
will comments on this approximation in the next section. 

An important feature exposed by (14) is that the cos 2(p an- 
gular dependence contributed via the Z boson has an opposite 
sign compared with the one contributed via a virtual photon. 
This can been seen from the fact that el > 0, while 



{Sv^ - Sl^) < for all flavors. 



(28) 



III. NUMERICAL RESULTS 

We will show that RHIC [56] at BNL is ideal for mea- 
suring the sign-reversal of v**'. The pp Drell-Yan process 
in the RHIC-spin program is mainly dedicated to polarized 
proton beams, although a spin-averaged measurement is still 
allowed. At RHIC, with the highest center of mass energy 
yfs = 500GeV, the experiments can probe the region x ~ 0.2 
(corresponding to the mid-rapidty of the vector boson) as the 
dilepton mass is near the Z-pole. We will present the predic- 
tion for 2v**', since it is equal to 2v + A- 1, and the later is 
the observable that can be directly measured by experiments. 

For the Boer-Mulders functions needed in the calculation, 
we adopt the parametrization [57] extracted from the unpolar- 
ized pd and pp Drell-Yan data at < (77- < 2 GeV measured 
by E866/NuSea Collaboration at FNAL. The E866/NuSea ex- 
periment covers the region 4.5 < Q < 15 GeV (exclud- 
ing the T resonance region) with a center of mass energy 
yfs - 38.7 GeV. Therefore the x region near the Z-pole at 
RHIC is similar to that of the E866/NuSea experiment. Also 
RHIC is a proton-proton collider where the dilepton is pro- 
duced by the annihilation of valence and sea quark from each 
hadron, just like the case of the E866/NuSea experiments. The 
kinematical cuts applied in the calculations are: 



-1.5<y<1.5, 0<?7-<2GeV, 



(30) 



' In principle there can be also a sin 2(p azimuthal dependence. However, it 
is l/Q suppressed compared to the cos 20 dependence and can be ignored 
here. 



where y = 5 ln(xi /x2) is the rapidity of the y* or Z boson. The 
reasons to choose a low pt cut are twofold. First at low qT the 
corrections from QCD are small, therefore the approximation 
in (27) is valid. The other is that in at this special low qT 
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FIG. 2: (a): The Q-dependent cos 20 asymmetry in pp — > -y*/Z + X — > {*{* — X process at RHIC. (b): The ^r-dependent cos 20 asymmetry in 
the unpolarized pp ^> y * /Z + X ^> y' /Z + X process at RHIC for diiTerent Q values 
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FIG. 3: (a): The 2-dependent cos 20 asymmetry in pp — > y* /Z + X — » £*£ + X process at RHIC. (b): The gj-dependent cos 20 asymmetry 
in the unpolarized pp -^ y'/Z + X -^ l*(^X process at RHIC for different Q values 



region the intrinsic transverse momentum of partons is more 
relevant, thus the tree-level result for v**^ is justified. 

In Fig. 2a we plot the prediction for v**^ (scaled with a fac- 
tor of 2, in order to correspond to the size of the Lam-Tung 
relation violation) as a function of dilepton mass Q at RHIC 
for ^fs - 500 GeV. The result clearly shows that the value 
of 2v**^ is positive as g < 60 GeV, while it reverses sign to 
be negative as 2 > 75 GeV, and reaches -5% at the Z-pole. 
In Fig 2b we plot the prediction for 2v^*^ as a function of 
transverse momentum of the dilepton qr, for Q - 50 GeV 
and Q - Mz- Again it shows a sign reversal of 2v**^ at 
two different Q value. The magnitude of the 2v**^ peak at 
Qt - IGeV falls at higher qr- The qj shape of the asym- 
metries indicate that the intrinsic transverse momentum of the 
parton plays a significant role at low qr, and can give a sub- 
stantial contribution. In Fig. 3a we plot the y-dependent re- 
sults of 2v*'^ for e = 50 GeV and Q = Mz, and in Fig. 3b 
we plot the y-dependent results of 2v^^ for 30 < 2 < 45 GeV 



and 75 < 2 < 110 GeV. An observation from these two 
figures is that the magnitude of 2v**^ increases as the rapidity 
increases. Our theoretical predictions suggests that the magni- 
tude of 2v-i-/l- 1 is sizable near the Z-pole and is measurable at 
RHIC. Therefore an accurate measurement on 2v + /I - 1 both 
at Z-pole and at a lower dilepton mass region can serve as a 
test of the chiral-odd property of quarks inside an unpolarized 
nucleon. 

Several points need to be addressed here. First, as the evolu- 
tion [58] of the TMD Boer-Mulders function still remains un- 
clear, we assume that the scale dependences of hf{x,pj) and 
the spin averaged distribution function /i (x, p^) are the same 
in calculating v**', which was also adopted before in the ex- 
traction of the Boer-Mulders function [57]. We would like to 
admit that the evolution of hf(x, pj) can be more complicated 
than that of /i (x, p^), as the former one is chiral-odd, while 
the later one is chiral-even. Further more, recent quantitative 



calculation [59] on the spin-independent processes demon- 
strate that the CoUins-Soper evolution of TMD distributions 
may be significant. Since v**^ is approximately the ratio be- 
tween hf and /i , the evolution eff'ect can only influence our 
results quantitatively at most, but not qualitatively. As our 
main purpose is to reveal the sign change of v^*' between 
the low Q region and Z mass region, our assumption on the 
scale dependence can be viewed as a reasonable choice. Sec- 
ondly, in Eqs. (10,11) we employ a tree-level expression of 
the TMD factorization formula, in which the soft factor has 
not been considered. Resummation effects [60-62] of soft 
gluon radiation in the TMD factorization will lead to a Sud- 
dakov factor that alters the tree-level result. The study in 
Ref. [61] shows that TMD azimuthal spin asymmetries are 
suppressed by this Suddakov factor in the region where qj 
is much larger than the intrinsic transverse momentum of the 
parton, but still much smaller than Q. In our calculation we 
restrict the cut on the transverse momentum of the dilepton as 
< qr < 2 GeV, where the intrinsic transverse momentum 
of partons plays a significant role, and we assume that the tree 
level approximation still holds, to avoid Suddakov suppres- 
sion. The Suddakov effect is certainly important for azimuthal 
observables at higher qj (but still much small than Q) and 
should be considered. Based on the uncertainties discussed 
above, our result can be viewed as an estimate. Nevertheless, 
our study provides a useful understanding of the cos 20 az- 
imuthal asymmetry in the Z mass region from the tree-level 
calculation. Very recently, new theoretical analysis given in 
Ref. [63] shows that, soft factors appearing beyond tree level 
cancel out of the weighted azimuthal asymmetry by employ- 
ing Bessel functions. We expect that this newer approach can 



be applied to study the weighted cos 2(f> azimuthal asymmetry 
at the Z-pole to provide a rigorous test on the sign change of 



IV. CONCLUSION 

We study the cos 20 angular dependence of dilepton pro- 
duction at the Z-pole. We show that, due to the chiral-odd na- 
ture of Boer-Mulders function, the Boer-Mulders effect will 
cause a sign change of 2v + /I - 1 (or equivalently, 2v**') from 
the low Q region to the Z mass region. This is a special signa- 
ture of the Boer-Mulders effect as a source of violation of the 
Lam-Tung relation. Therefore, the experimental detection of 
the sign reversal of 2v + A - 1 will serve as a clear evidence 
of the existence of chiral-odd quarks inside an unpolarized nu- 
cleon. Using a recent extracted set of Boer-Mulders functions, 
we predict the sign and magnitude of v**' near the Z-pole at 
RHIC with V^ = 500 GeV, based on TMD factorization. Our 
analysis at the low Q and Z mass regions shows that this sign 
reversal can be detected if the Boer-Mulders effect indeed is 
a source of Lam-Tung relation violation. A dedicated mea- 
surement on both of the size and the sign of the Lam-Tang 
violation at RHIC in a wide Q range therefore provides valu- 
able hint on the dynamics of the cos 20 azimuthal dependence 
in dilepton production. 

Acknowledgements This work is supported by FONDE- 
CYT (Chile) Project Nos. 11090085, No. 1100715, by Project 
Basal FB0821, and by NSFC (China) Project No. 1 1005018. 



[1] For a recent review on the related topics, see V. Barone, 
F. Bradamante, and A. Martin, Prog. Part. Nucl. Phys. 65, 267, 
(2010). 
[2] R. J. Cakes, Nuovo Cimento 44A, 440 (1966). 
[3] C. S. Lam and W. K. Tung, Phys. Rev. D 18, 2447 (1978). 
[4] C. S. Lam and W. K. Tung, Phys. Lett. B 80, 228 (1979). 
[5] C. S. Lam and W. K. Tung, Phys. Rev. D 21, 2712 (1980). 
[6] E. Mirkes and J. Ohnemus, Phys. Rev. D 50, 5692 (1994); Phys. 

Rev. D 51, 4891 (1995). 
[7] D. Boer and W. Vogelsang, Phys. Rev. D 74, 014004 (2006). 
[8] E. L. Berger, J. Qiu, and R. A. Rodriguez-Pedraza, Phys. Lett. 

B 656, 74 (2007); Phys. Rev. D 76, 074006 (2007). 
[9] S. Falciano, et al, NAIO Collaboration, Z. Phys. C 31, 513 
(1986); M. Guanziroli, et al, NAIO Collaboration, Z. Phys. C 
37, 545 (1988). 
[10] J.S. Conway, et al, Phys. Rev. D 39, 92 (1989). 
[11] L.Y. Zhu, et. al, (FNAL-E866/NuSea Collaboration), Phys. 

Rev. Lett. 99, 082301 (2007). 
[12] L. Y. Zhu et al (FNAL E866/NuSea Collaboration), Phys. Rev. 

Lett. 102, 182001 (2009). 
[13] T. Aaltonen et al (CDF Collaboration), Phys. Rev. Lett. 106, 

241801 (2011). 
[14] A. Brandenburg, O. Nachtmann, and E. Mirkes, Z. Phys. C 60, 

697 (1993). 
[15] D. Boer, A. Brandenburg, O. Nachtmann, and A. Utermann, 
Eur. Phys. J. C 40, 55 (2005). 



[16] A. Brandenburg, S.J. Brodsky, V.V. Khoze, and D. Miiller, 

Phys. Rev. Lett. 73, 939 (1994). 
[17] K.J. Eskola, P. Hoyer, M. Vanttinen, and R. Vogt, Phys. Lett. 

B333, 526 (1994). 
[18] D. Boer, Phys. Rev. D 60, 014012 (1999). 
[19] D. Boer and RJ. Mulders, Phys. Rev. D 57, 5780 (1998). 
[20] Z. Lu and B.-Q. Ma, Phys. Rev. D 70, 094044 (2004). 
[21] Z. Lu and B. Q. Ma, Phys. Lett. B 615, 200 (2005). 
[22] A. Bianconi and M. Radici, Phys. Rev. D 72, 074013 (2005). 
[23] A. N. Sissakian, O. Yu. Shevchenko, A.P. Nagaytsev, and O.N. 

Ivanov, Phys. Rev. D 72, 054027 (2005). 
[24] A. Sissakian, O. Shevchenko, A. Nagaytsev, O. Denisov, and 

O. Ivanov, Eur Phys. J .C 46, 147 (2006). 
[25] Z. Lu, B. Q. Ma and I. Schmidt, Phys. Lett. B 639, 494 (2006). 
[26] V. Barone, Z. Lu and B. Q. Ma, Eur. Phys. J. C 49, 967 (2007). 
[27] Z. Lu, B. Q. Ma and I. Schmidt, Phys. Rev. D 75, 014026 

(2007). 
[28] Z. Lu, B.-Q. Ma, and I. Schmidt, Phys. Rev. D 75, 014026 

(2007). 
[29] L.R Gamberg and G.R. Goldstein, Phys. Lett. B 650, 362 

(2007). 
[30] B. Zhang, Z. Lu, B.-Q. Ma and I. Schmidt, Phys. Rev. D 77, 

054011 (2008); ibid. 78, 034035 (2008). 
[31] V. Barone, S. Melis and A. Prokudin, Phys. Rev D 82, 114025 

(2010) 
[32] J.C. Colhns and D.E. Soper, Phys. Rev. D 16, 2219 (1977). 



[33] S.J. Brodsky, D.S. Hwang, and I. Schmidt, Phys. Lett. B 530, 

99 (2002). 
[34] J.C. Collins, Phys. Lett. B 536, 43 (2002). 
[35] D. Boer, S.J. Brodsky, and D.S. Hwang, Phys. Rev. D 67, 

054003 (2003). 
[36] A.V. Belitsky, X. Ji, and F. Yuan, Nucl. Phys. B656, 165 (2003). 
[37] D. Boer, PJ. Mulders, and F. Pijlman, Nucl. Phys. B667, 201 

(2003). 
[38] G. R. Goldstein and L. P. Gamberg, in Proceedings of the 

31st International Conference on High Energy Physics (ICHEP 

2002), Amsterdam, 2002 (Elsevier Science, Amsterdam, 2002), 

pp. 452C454. 
[39] PV. Pobyhtsa, hep-ph/0301236. 
[40] F Yuan, Phys. Lett. B 575, 45 (2003). 
[41] A. Bacchetta, A. Schafer, and J.-J. Yang, Phys. Lett. B 578, 109 

(2004). 
[42] M. Burkardt and B. Hannafious, Phys. Lett. B 658, 130 (2008). 
[43] L. P Gamberg, G. R. Goldstein and M. Schlegel, Phys. Rev. D 

77, 094016 (2008). 
[44] S. Meissner, A. Metz, and K. Goeke, Phys. Rev. D 76, 034002 

(2007). 
[45] A. Bacchetta, F Conti and M. Radici, Phys. Rev. D 78, 074010 

(2008). 
[46] M. Wakamatsu, Phys. Rev. D 79, 094028 (2009). 
[47] A. Courtoy, S. Scopetta and V. Vento, Phys. Rev. D 80, 074032 

(2009). 



[48] A. Bacchetta, M. Radici, F Conti and M. Guagnelli, Eur. Phys. 
J. A 45, 373 (2010). 

[49] B. Pasquini and F. Yuan, Phys. Rev. D 81, 1 14013 (2010). 

[50] B. Pasquini and P Schweitzer, Phys. Rev. D 83, 1 14044 (201 1). 

[51] M. Gockeler et al (QCDSF Collaboration and UKQCD Col- 
laboration), Phys. Rev. Lett. 98, 222001 (2007). 

[52] S. Meissner, A. Metz, M. Schlegel and K. Goeke, J. High En- 
ergy Phys. 08 (2008) 038. 

[53] L. Gamberg and M. Schlegel, Phys. Lett. B 685, 95, 2010. 

[54] D. Brommel et al. (QCDSF Collaboration and UKQCD Collab- 
oration), Phys. Rev. Lett. 101, 122001 (2008). 

[55] S. Arnold, A. Metz, and M. Schlegel, Phys. Rev. D 79, 034005 
(2009). 

[56] G. Bunce, N. Saito, J. Soffer and W. Vogelsang, Ann. Rev. Nucl. 
Part. Sci. 50, 525 (2000). 

[57] Z. Lu and L Schmidt, Phys. Rev. D 81, 034023 (2010). 

[58] A. A. Henneman, D. Boer, and P J. Mulders, Nucl. Phys. B620, 
331-350(2002). 

[59] S. M. Aybat, T. C. Rogers, Phys. Rev. D 83, 1 14042 (201 1). 

[60] J. C. Collins, Nucl. Phys. B396, 161 (1993). 

[61] D. Boer, Nucl. Phys. B603, 195 (2001). 

[62] Z. B. Kang, B. W. Xiao and F. Yuan, arXivil 106.0266. 

[63] D. Boer, L. Gamberg, B. Musch, A. Prokudin, 
arXiv: 1107.5294. 



